fun = @(x,y) (y.^2.*exp(-y.^2)+x.^4.*exp(-x.^2))./x.*exp(-x.^2);
q = integral2(fun,2,4,-1,1)